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Abstract. We discuss observations of the linearly polar- 
ized component of the diffuse galactic radio background. 
These observations, with an angular resolution of 4', were 
made with the Westerbork Synthesis Radio Telescope 
(WSRT) in 5 frequency bands in the range 341 - 375 
MHz. The linearly polarized intensity P (with polarized 
brightness temperature going up to 10 K) shows a 'cloudy' 
structure, with characteristic scales of 15 - 30', which con- 
tains relatively long, but very narrow 'canals' (essentially 
unresolved) in which P is only a small fraction of that in 
the neighbouring beams. 

These 'canals' are generally seen in more than one fre- 
quency band, although their appearance changes between 
bands. They are probably due to depolarization within 
the synthesized beam, because the change in polarization 
angle A0 po i across the deepest 'canals' is in general close 
to 90° (or 270° etc.). These very abrupt changes in 4> po i, 
which are seen only across the 'canals', seem to be ac- 
companied by abrupt changes in the Rotation Measure 
(RM), which may have the right magnitude to create the 
difference of close to 90° in 4> po \, and thereby the 'canals'. 

The structure in the polarization maps is most likely 
due to Faraday rotation modulation of the probably 
smooth polarized radiation emitted in the halo of our 
Galaxy by the fairly local (< 500 pc) ISM. Therefore, 
the abrupt changes of RM across the 'canals' provide evi- 
dence for very thin (^ 1 pc), and relatively long transition 
regions in the ISM, across which the RM changes by as 
much as 100%. Such drastic RM changes may well be due 
primarily to abrupt changes in the magnetic field. 
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1. Introduction 

Wieringa et al. (1993) were the first to note structure on 
arcminute scales in the linearly polarized component of 
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the galactic radio background at 325 MHz, observed with 
the WSRT. The small-scale structure in the maps of polar- 
ized intensity P, (with polarized brightness temperatures 
Tb iPO i of up to 10 K) does NOT have a counterpart in 
total intensity, or Stokes /, down to very low limits. Be- 
cause the total Stokes I of the galactic radio background 
has an estimated Tb. po i of the order of 30 - 50 K at 325 
MHz, which must be very smooth and therefore filtered 
out completely in the WSRT measurements, the appar- 
ent polarization percentage of the small-scale features can 
become very much larger than 100%. 

The absence of corresponding small-scale structure 
in Stokes I led Wieringa et al. (ibid.) to propose that 
the small-scale structure in polarized intensity P is due 
to Faraday rotation modulation. In this picture, syn- 
chrotron radiation generated in the Galactic halo reaches 
us through a magneto-ionic screen, viz. the warm rel- 
atively nearby ISM. Structure in the electron density 
and/or magnetic field in the ISM causes spatial variations 
in the Rotation Measure (RM) of the screen. Hence, the 
angle of linear polarization of the synchrotron emission 
from the halo is rotated by different amounts along dif- 
ferent lines of sight. Even if the polarized emission in the 
halo were totally smooth, in intensity as well as angle, the 
screen would produce structure in Stokes Q and U . 

Small-scale structure in the polarized galactic radio 
background has recently been observed also at other fre- 
quencies. At 1420 MHz, Gray et al. (1998, 1999) used the 
DRAO synthesis telescope to study the phenomenon at 
1' resolution. Uyaniker et al. (1999) used the Effelsberg 
telescope at 1.4 GHz, to map the polarized emission at 9' 
resolution over about 1100 D . Duncan et al. (1998) discuss 
radio polarization data at 1.4, 2.4 and 4.8 GHz with the 
Parkes radio telescope and the VLA, at resp. 5', 10' and 
15' resolution. All these observations support the inter- 
pretation in terms of modulation of emission originating 
at larger distances, by a relatively nearby Faraday screen. 

The distributions of polarized intensity and angle may 
therefore be used to study the structure of the Faraday 
screen. In particular, polarization observations give infor- 
mation about the electron density, n c , and the component 
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of the magnetic field parallel to the line of sight, B\\ , in 
the ISM on scales down to less than ~ 0.5 pc (< 4' at 
an assumed distance of ~ 500 pc). The diffuse nature of 
the polarized radio background allows (almost) complete 
spatial mapping of RMs over large areas, provided one 
has observations at several frequencies. This gives a large 
advantage over RM determinations through individual ob- 
jects, like pulsars or extra-galactic radio sources. 

2. Distribution of polarized intensity 

In Fig. |l| we show a gray-scale representation of the polar- 
ized intensity in a 5 MHz wide frequency band centered 
at 349 MHz. The map shows a region of 6.4° x 9° centered 
at a = 6' l 10 m ,<5 = 53°(£ = 161°, b = 16°) at an angu- 
lar resolution of about 4'. It is one of 8 frequency bands 
observed simultaneously. Three of those have strong inter- 
ference, but we obtained good data at 341, 349, 355, 360 
and 375 MHz. All 5 maps were made combining mosaics of 
7x5 pointing centres. This yields constant sensitivity over 
a large area (see e.g. Rengelink et al. 1997). The observa- 
tions were made with the WSRT in January and February 
1996, largely at night, and ionospheric Faraday rotation 
was therefore well-behaved. No corrections were applied. 

The region in Fig. |l| is rather special because Tb iPO i 
goes up to 10 K, and because it contains large, almost 
linear structures in P. Our attention was drawn to this 
field by the panoramic view of galactic polarization pro- 
duced in the WENSS survey (de Bruyn & Katgert 2000). 
However, this field is not unique, and there are other re- 
gions with similarly high Tb jPO i. Over a very large fraction 
of the map the P-signal is quite significant, with a noise 
er Tb w 0.5 K. With S/N-ratios of generally more than 3 
and going up to 30, polarization angles are well-defined. 
Note that in this region, the upper limit to structure in 
Stokes I (total intensity) on small scales (^ 30') is about 
1 K, or less than 2% of the total /. 

There appear to be at least two distinct components 
in the polarized intensity distribution. The first one is a 
fairly smooth, 'cloudy' component, pervading the entire 
map, with intensity variations on typical scales of (sev- 
eral) tens of arcminutes. In addition, there are conspicu- 
ous, very narrow and often quite long and wiggly struc- 
tures, which we will refer to as 'canals', in which the po- 
larized intensity is considerably lower than in the imme- 
diate surroundings. In this Letter we focus on the nature 
and implications of the narrow 'canals'; we will discuss 
the 'cloudy' component in more detail in another paper 
(Haverkorn et al. 2000). 

3. The nature of the 'canals' in polarized intensity 

The strong and abrupt decrease of polarized intensity in 
the 'canals' suggests that depolarization is responsible. 
There are several mechanisms that can produce depolar- 
ization, but the only plausible type in this case is beam 
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Fig. 1. Linearly polarized intensity P at 349 MHz in a 
6.4° x 9° field centered at I = 161°, 6=16°. The resolution 
is ~ 4', the maximum brightness temperature is ~ 10 K. 
The generally 'cloudy' distribution contains long narrow 
'canals' of low P. The white box shows the area displayed 
in Fig. |. 

depolarization. This occurs when the polarization angle 
varies significantly within a beam. Complete depolariza- 
tion requires that for each line of sight there is a 'com- 
panion' line of sight within the same beam that has the 
same polarized intensity but for which the polarization 
angle differs by 90° . Below we will show that our observa- 
tions indicate that the polarization angle indeed changes 
by large amounts across low polarized intensity 'canals', 
and close to 90° across the 'canals' of lowest P. 

Depolarization can also be caused by 'differential Fara- 
day rotation'. This happens when along a line of sight 
emitting and (Faraday) rotating plasmas coexist (e.g. 
Burn 1966; Sokoloff et al. 1998). However, the absence 
of correlated structure in Stokes / and the high degree 
of polarization suggest that this is not a dominating ef- 
fect. Significant bandwidth depolarization, which occurs 
when the polarization angle is rotated by greatly different 
amounts in different parts of a frequency band could only 
play a role (given our 5 MHz bandwidth) if the RM were 
of order 80 rad m~ 2 , which is not the case in this region 
near the galactic anti-centre (see below). 

In Fig. H we show the polarization vectors around a few 
of the deepest 'canals', superimposed on gray-scale plots 
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Fig. 2. Polarized intensity P at 349 MHz (left) and 
360 MHz (right) of the area inside the box in Fig. ||. Po- 
larization angles and intensities are indicated by the vec- 
tors, which are sampled at locations 4' apart (independent 
beams). Note A</> po i w 90° across low-P 'canals'. 

of P, in two frequency bands. The area shown is indicated 
in Fig. |l|. The polarization vectors on either side of the 
'canals' are quite close to perpendicular, demonstrating 
that the 'canals' are produced by beam depolarization. 
This perpendicularity applies to all 'canals', irrespective 
of frequency band and is very convincing, especially be- 
cause everywhere else the polarization vectors vary quite 
smoothly (if significantly!). 

Beam depolarization creates 'canals' that are one beam 
wide, which is exactly what we observe. This implies that 
the 90° 'jump' must occur on angular scales smaller than 
the beamwidth. At ~ 2' resolution (about twice that in 
Fig. ||), the 'canals' indeed seem unresolved, but the de- 
crease in S/N-ratio precludes conclusions on even smaller 
scales (the original data have 0.8' resolution). 

Additional evidence that the 'canals' are due to beam 
depolarization is statistical. We defined 'canal-like' points 
from the observed values of P, as follows. For each point 
in the mosaic we compared the observed value of P with 
the P-values in pairs of two diametrically opposed neigh- 
bouring (adjacent) points. If the value of P in the central 
point was less than a certain small fraction of the values in 
both comparison points, the point was defined 'canal-like'. 
This definition mimics the visual detection 'algorithm'. 

In the top panel of Fig. || we show the distribution 
of the difference between the 4> po \s in the two adjacent 
points that define the 'canal-like' points, for a P-threshold 
of 30%. The A0 po i-distribution peaks at 90°, fully consis- 
tent with the beam depolarization hypothesis. This con- 
clusion is reinforced by a comparison with the distribution 
of A</>p i (again for diametrically opposed adjacent neigh- 
bours) of all points for which P is between 1.0 and 2.0 
times larger than both P-values in the two diametrically 
opposed neighbouring points, which is shown in the bot- 
tom panel of the same figure. 

Similar 'canals' were noted by Uyaniker et al. (1999) 
and Duncan et al. (1998), who also invoked beam depolar- 
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Fig. 3. Top panel: Polarization angle difference A0 po i be- 
tween two points on opposite sides of a 'canal-like' point 
(see text for definition) . A clear preference for A(/> po i ~ 90° 
across 'canals' is visible. 

Bottom panel: A</> po i between two points on opposite sides 
of non-'canal-like' points (see text for definition). For the 
non-'canal-like' points , (A</> po i) as 0°, rather than 90°. 

ization. Yet, Fig. [| is the first quantitative proof for this 
explanation. 

4. The cause of the 'jumps' in polarization angle 

Two processes can cause jumps in polarization angle 4> po \ 
across the 'canals': a sudden change in RM across the 
'canals', and a jump in intrinsic po i of the emission in- 
cident on the Faraday screen. A large change in intrinsic 
</> po i implies a change in magnetic field direction and is 
therefore quite difficult to understand in view of the ab- 
sence of structure in total intensity / at the more than 2% 
level (see Sect. ||). On the other hand, variations in the RM 
of the Faraday screen would seem to be quite natural, if 
not unavoidable. 

Discontinuities in RM must play an important role in 
producing the 'canals', because the 'canals', although sim- 
ilar in adjacent frequency bands, generally do not occur in 
all bands, and certainly are not identical in the different 
bands (see Fig. ||). This indicates that the jumps in 4> po \ 
are mainly due to changes in RM. However, the question 
is if the jumps in </> po i are indeed accompanied by jumps 
in RM of the right magnitude so that A<p po \ — 90° is pro- 
duced at the frequency where the 'canal' is best visible. 

In principle, the determination of RM only involves 
a simple linear fit of the polarization angles in the five 
frequency bands (at 341, 349, 355, 360 and 375 MHz) vs. 
A 2 , but in practice several complications may arise. First, 
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Fig. 4. Polarized intensity P at 349 MHz in the area in- 
dicated by the box in Fig. g. Overlaid are small plots of 
^poi(A) vs. A 2 for independent points, with a linear fit 
through the data. RM's range from ~ -7 rad m -2 in the 
upper left corner to -1.5 rad m~ 2 in the 'island' between 
the two 'canals' (below centre). 



the observed values of </> po i may be biased due to imaging 
effects (like off-sets) in the Stokes Q- and /7-maps from 
which 4> po \ is derived (cf. Wieringa et al. 1993). Our data 
indicate that, in the maps of this region of sky, such off- 
sets are quite small, so that the bias in the 4> po \ values 
is small. Second, it is not obvious that the assumption of 
pure Faraday rotation (0(A) oc A 2 ) is supported by the 
data (see Haverkorn et al. 2000). 

In Fig. H we show an array of plots of po i(A) vs. A 2 
for independent beams in the small region (indicated in 
Fig. ||) that contains two clear 'canals'. As can be seen, a 
direct determination of ARM across the 'canals' is not at 
all trivial. Without knowing the position of the 'canals', 
one probably would have some trouble to find the 'canals' 
from discontinuities in RM distribution alone, due to the 
uncertainties in the RM-cstimates, which sometimes are 
considerable. On the other hand, if one knows where the 
canals are one can identify some related 'jumps' in RM. 

From the present data, it seems quite likely that the 
'canals' are primarily due to quite abrupt and relatively 
large changes of RM, with ARM/RM ranging from ~ 0.3 
to more than 1 (at least in this region of sky). Note that 
in this region the RMs are in the range from -10 to +10 
rad m~ 2 (also confirmed by several polarized extragalac- 
tic radio sources in these same observations). However, a 



more robust conclusion about the relation between A0 po i 
and ARM requires a detailed analysis of more, and more 
sensitive data, and a careful error analysis. 

5. Implications for the structure of the local ISM 

Because we have not yet reached a quantitative conclusion 
about the suspected correlation between A(/) po i and ARM, 
it is not possible to give a full discussion of the implica- 
tions that these polarization data have for the small-scale 
structure of the warm ISM. However, the data discussed 
here show the great promise that high-resolution, multi- 
band polarization data hold for the study of the ISM, es- 
pecially on small scales where pulsars and extragalactic 
radio sources cannot give much information. 

Fortunately, more and more sensitive radio polariza- 
tion data (in different regions of sky) are forthcoming. In 
addition, information must be obtained about the elec- 
tron density in the warm ISM on the relevant scales (e.g. 
through Ha measurements), as well as on the other com- 
ponents in the ISM (like e.g. the HI). 

While wc fully realize the preliminary nature of the 
conclusions presented, we feel justified to speculate some- 
what on the possible implications of the 'canals'. Structure 
in RM reflects structure in B\\ and/or n c in the ISM. How- 
ever, as the RM is an integral over the entire line of sight, 
the large ARM/RM values that are implied by our ob- 
servations may give a very specific message. In particular, 
we consider it unlikely that the large ARM/RM values are 
produced mainly by variations in electron density. Instead, 
they may be indicating a turbulent ISM with varying (re- 
versing) magnetic field structures, as modeled in recent 
MHD simulations (see e.g. Mac Low & Ossenkopf 2000; 
Vazquez-Semadeni & Passot 1999). 

Acknowledgements. The Westerbork Synthesis Radio Tele- 
scope is operated by the Netherlands Foundation for Research 
in Astronomy (NFRA) with financial support from the Nether- 
lands Organization for scientific research (NWO). This work is 
supported by NWO grant 614-21-006. 

References 

Burn B.J., 1966, MNRAS 133, 67 
de Bruyn A.G., Katgert P., 2000, in preparation 
Duncan A.R., Haynes R.F., Reich W. et al., 1998, MNRAS 
299, 942 

Gray A.D., Landecker T.L., Dewdney P.E., Taylor A.R., 1998, 
Nat 393, 660 

Gray A.D., Landecker T.L., D ewdney P.E. et al., 1999, ApJ 
514, 221, |astro-ph/9812164| 



Haverkorn M., Katgert P., de Bruyn A.G., 2000, in prepa ration 
Mac Low M Ossenkopf V., 2000, A&A 353, 339 
| ph/9906334] 



astro- 



Rengelink, R.B., Tang Y., de Bruyn A.G. et al., 1997, A&AS 
124, 259 

Sokoloff D.D., Bykov A.A., Shukurov A., et al., 1998, MNRAS 
299, 189 



M. Haverkorn, P. Katgert, A.G. de Bruyn: Structure in the local Galactic ISM on scales down to 1 pc 



5 



Uyaniker B., Fiirst E., Reich W., Reich P ., Wielebinski, R., 
1999, A&AS 138, 31, [astro-ph/990502; 

Vazquez- Semadeni E., Passot T., 1999, in Interstellar Turbu- 
lence, Franco J., Carraminana, A., Cambridge Univ. Press, 
p. 223 

Wieringa M.H., de Bruyn A.G., Jansen D., Brouw W.N., Kat- 
gert P., 1993, A&A 268, 215 



